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Abstract 

The Baldwin-Lomax algebraic turbulence model was 
modified for hypersonic flow conditions. Two coefficients 
in the outer-layer eddy-viscosity model were determined 
as functions of Mach number and temperature ratio. By 
matching the solutions from the Baldwin-Lomax model to 
those from the Cebeci-Smith model for a flat plate at hy- 
personic speed, the new values of the coefficients were ob- 
tained. The results show that the values of C ep and Ckicb 
are functions of both Mach number and wall temperature 
ratio. The C cp and Ckub variations with Mach number and 
wall temperature were used for the calculations of both a 
4° wedge flow at Mach 18 and an axisymmetric Mach 20 
nozzle flow. The Navier-Stokes equations with thin-layer 
approximation were solved for the above hypersonic flow 
conditions and the results were compared with existing 
experimental data. The agreement between the numerical 
solutions and the existing experimental data were good. 
The modified Baldwin-Lomax model thus is useful in the 
computations of hypersonic flows. 

Introduction 

With recent interest in hypersonic flight, the flow phe- 
nomena at high Mach number are gaining the attentions 
of both, experimental and computational fluid dynamics 
communities. There has been much progress in turbulence 
modeling for incompressible flows, but little progress in 
high speed flow. For hypersonic turbulent boundary layer 
calculations, most investigators have used algebraic turbu- 
lence models for their simplicity, and a popular model is 
the Baldwin-Lomax model 1 which was patterned after the 
Cebeci-Smith model 2 . It uses the same inner- layer eddy- 
viscosity formula as the Cebeci-Smith model, but it does 
not require the finding of the edge of the boundary layer 
for the outer-layer eddy-viscosity calculation. 

There are two coefficients in the outer-layer eddy- 
viscosity formula of the Baldwin-Lomax turbulence model 
for attached wall boundary layer. Their values were ob- 
tained by matching the outer-layer eddy-viscosity with 
that from the Cebeci-Smith model at transonic 


flow conditions 1 . Their suggested values at these condi- 
tions are 1.6 for C cp and 0.3 for Cktcb- Several authors 
have suggested different values for the coefficients, depend- 
ing on Mach number and pressure gradient( 3,4,6,e ). As 
free-stream Mach number increases, the sublayer thickness 
increases and this growth of sublayer becomes significant 
at hypersonic Mach numbers. To obtain the same perfor- 
mance of the Baldwin-Lomax model as that of the Cebeci- 
Smith model, the two coefficients should be adjusted with 
Mach number and wall temperature ratio. 

Previous investigators have examined the variation of 
C cp and Ckicb with pressure gradient, Mach number, and 
surface blowing. Table 1 lists variations of C cp from 0.85 
to 1.82 and variations of Ckicb from .3 to .65. Granville 
determined the pressure gradient effects on both of the co- 
efficients for incompressible flow. York and Knight found 
that both the coefficients depend on local Mach number, 
local skin-friction, and the pressure parameter. York and 
Knight’s investigation included the flows up to Mach 5. 
The present authors did not find in the literature any eval- 
uations of the coefficients in the Baldwin-Lomax model at 
hypersonic speeds. The above studies suggest that C cp and 
CkUb should be adjusted with the local flowfield. However, 
at the present time it is typical to use fixed coefficients in 
the Baldwin-Lomax turbulence model. It is an objective of 
this study to determine the variation of these coefficients 
with Mach number and wall temperature ratio. 

To study the variation of the coefficients in the 
Baldwin- Lomax model, the boundary layer over a flat 
plate was solved with the Baldwin-Lomax model and 
the results were compared with those from the Cebeci- 
Smith model. By matching the skin friction coefficients 
and boundary layer thicknesses from the Baldwin- Lomax 
model to those from the Cebeci-Smith model, new coef- 
ficients were obtained at different Mach numbers. These 
new coefficients were used for the calculations of the flows 
over a wedge and in a nozzle at hypersonic Mach numbers, 
and the results were compared with experimental data. 


Turbulence models 
Cebeci-Smith model 

The Cebeci-Smith model can be written as follows: 
Ht = min(nu,Vto) (l) 
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where 5,*,,,. is the incompressible form of boundary layer 
displacement thickness, A + is the van Drlest’s damping 
factor(= 26) and ki is constant{= .0168). 


Baldwin-Lomax model 

The Baldwin-Lomax model can be written as follows: 

Ht = min{nti,nto) ( 6 ) 

Pu = pi 1 M (7) 

IHo = C cp p F wa kt Fkub (8) 
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where C cp is 1.6, Ckitb is .3, C w k is 0.26, and u is the 
vorticity. The values of y max and F max are obtained from 
the following function: 
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and the value of Uhj is the difference between the maxi- 
mum and minimum total velocities In the profile. 

As seen from the above formulas, the Baldwin-Lomax 
model is nearly same as the Cebeci-Smith model for the 
inner-layer. For outer-layer formulas, the following rela- 
tions can be obtained by matching the above two models: 
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Ckub = 


Vmax 
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As seen from Eq. 12, C cp is a function of the incompress- 
ible form of the displacement thickness. As Mach number 
increases, the contribution to S' nc from the viscous sub- 
layer becomes significant. Unless y max and F max contain 
the effect of the growth of the viscous sublayer, C ep is func- 
tions of both Mach number and wall temperature ratio. 


To study the variations of C cp and Ckub with Mach 
number and wall temperature ratio, a turbulent boundary 
layer over a flat plate was solved by the STAN5 program 7 . 
By matching boundary thickness and skin friction coeffi- 
cient from the Baldwin-Lomax model to those from the 
Cebeci-Smith model, C cp and Ckitb variations with Mach 
number and wall temperature ratio were obtained. For 
simplicity, the Klebanoff intermittency factor was assumed 
to be unity because main interest is in hypersonic flow 8 . 

Figure 1 shows the variation of C cp with Mach number 
for the following conditions. 

Flat plate 

Adiabatic wall, helium 
T t = 7°K 

R ex = 9.4 x 10 7 

For the calculation of laminar viscosity of helium, the fol- 
lowing viscosity laws were used: 

©(*■<*’*) <h) 

n = 5.0277 x 10~ 7 T 0 ' 647 (15) : 

and the turbulent Prandtl number was assumed to be con- 
stant, 0.9. 

As seen from Figure 1, C cp increases from 1.5 to 3.0 
as Mach number increases from 0 to 13. This increase in 
C cp with free-stream Mach number is mainly due to the 
growth of sublayer. Figure 2 shows the variation of Ckub 
with Mach number at the above condition. As seen from 
this Figure, Ckitb decreases with Mach number up to M=5 
and stays nearly constant. This variation shows that 
is nearly constant at hypersonic Mach number range. 

Figure. 3 and 4 show the variations of C cp and Ckub 
with Mach number for fully developed turbulent boundary 
layer over a flat plate for the following conditions: 

Flat plate 

Helium. 

Adiabatic wall and p- = .6 

To, = 500° K, Pc, = 270 atm. 

In this case, the free-stream conditions were obtained from 
isentropic relations and the Reynolds number range is from 
. 1.4 x 10 8 to 17 x 10 s . 

Figure 3 shows the variations of C cp with Mach num- 
ber. As seen from this Figure, C cp increases from about 
1,5 at Mach 1 to about 3.0 at Mach 20. The C ap values for 
the adiabatic conditions are greater than those for the cold 
wall conditions at Mach numbers higher than 10. This dif- 
ference is due to the fact that the viscous sublayer at the 
adiabatic wall conditions is thicker than that at the cold 
wall conditions. 

Figure 4 shows the variations of Ckitb with Mach num- 
ber. As seen from this Figure, Ckitb decreases from about 
0.6 at Mach 1 to about 0.37 for the adiabatic wall and 0.3 
for the cold wall at Mach 20. The Ckitb values for the adi- 
abatic wall conditions were nearly same as those for the 
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cold wall conditions up to Mach 8 and the values for 
adiabatic wall conditions are greater than those for cold 
wall conditions at very high Mach number. 

The above study indicates that C ep and C k i,b depend 
on Mach number and wall temperature ratio and vary from 
1.5 to 3.1 and 0.56 to 0.3, respectively. These variations of 
C cp and Ckub were used for the calculation of the hyper- 
sonic flows over a wedge and in a nozzle. 


Numerical Method and Grid 

An existing computer code, PARC2D, is used for the 
calculations for hypersonic Sows over a wedge and in a noz- 
zle. The code solves the two-dimensional Reynolds aver- 
aged Navier-Stokes equations with the Beam and Warming 
algorithm. The code was originally developed as AIR2D 
by Pulliam and Steger 9 , and it was called ARC2D after 
the Jameson type artificial dissipation 10 was added by 
Pulliam 11 . Later, Cooper adapted the ARC2D code for 
propulsion applications and named the code PARC2D 12 . 

Algebraically generated grids are used for the calcu- 
lations, The meshes are clustered near the wall with 100 
points vertically. The computational domain consists of 
uniform grid spacing up to x = lm and stretched grid 
spacing when x > lm in the streamwise direction, and the 
grid contains 150 x 100 points for the wedge flow calcu- 
lations. For the nozzle flow calculations, clustered grids 
in the streamwise direction are used near the throat. The 
grid size is 100 x 100 for the nozzle flow. 


Experimental Data 

To compare the numerical solutions with the experi- 
mental data, the following two cases were selected. 

1) Hypersonic flow over 4° wedge 13 . 

This experiment was conducted in the Mach 20 leg of the 
Langley high Reynolds number helium tunnel by Watson 
and the test conditions are as follows: 

Helium. 

Pt = 136 atm 
T, = 311 0 /f 
Mi = 18.05 
= .95 

= 4.36 x 10 7 /m 


T, 

T, 

EmI 

m 


Transition starts at x = .36 m and ends at x = ,89 m. For 
this wedge flow case, the shear stress at the wall, velocity 
and temperature profiles are known. The numerical so- 
lutions were compared with wall shear stress distribution 
and with velocity and temperature profiles at x = 2.159 
m. 

2) Axisymmetric contoured hypersonic nozzle flow 14 . 

For this case, the experimental data were obtained by 
Kemp and Owen in the Ames M-50 helium tunnel. The 
test condition is as follows: 

Helium. 

Pj = 270 aim. 


T t = 500 °K 
T w = 300 °K 
M tz it = 20 

For this case, the velocity and temperature profiles at 
x = 0.508 were compared with the numerical solutions. 


Modifications to Baldwln-Lomax Model 

As previously mentioned, the values of C cp and CkUb 
obtained for a flat plate were used for the calculations. 
Additionally, the following modifications were added to the 
Baldwin-Lomax model. 

1) Intermittency factor 

Experimental data show that the intermittency factor for 
hypersonic turbulent boundary layer is quite different from 
the incompressible one^ 8,13 ^. Instead of showing the vari- 
ation of Eq. 5, the intermittency factor stays nearly unity 
within the boundary layer and drops very sharply to zero 
at the edge of the boundary layer. To include the above 
intermittency variation within the boundary layer, the in- 
termittency factor is assumed to be unity up to 6 = 
and drop to zero at the edge of the boundary layer. 

2) Low Reynolds number effect 

To include the low Reynolds number effect on ki in the 
outer layer formula, Eq. 8, the following modification was 
used 15 ; 

, .0168(1 + 0.55) 

(i + n) 

n = .55(1 - exp(— 0.2432’ - 0.2982)) 

<«) 

where II is the strength of the wake and R e $ is the momen- 
tum thickness Reynolds number of the boundary layer. To 
calculate the momentum thickness of boundary layer, it 
was assumed that the point where y is is the edge of 
the boundary layer. 

3) Pressure gradient effect on A + 

To include the pressure gradient effect, the following modi- 
fied van Driest’s damping constant by Launder 16 was used: 


A + = 


26 


(£) 


1.7 


(15) 


where t is the local shear stress. 

Comparison with Experimental Data 

Calculations were carried out for the above hypersonic 
flow conditions with PARC2D code with thin-layer approx- 
imation using the Baldwin-Lomax turbulence model with 
the above modifications. 

1) Hypersonic wedge flow 

To calculate both laminar and turbulent boundary 
layers, the turbulent viscosity was turned on at the mid- 


3 



point between the beginning point and the end point of 
the experimentally determined transition (x = .61 m),e.g. 
instantaneous transition was assumed. 

Figure 5 shows the Mach number contour for the 4° 
wedge flow at free-stream Mach number of 18.05. The 
shock wave and boundary layer can be observed from this 
figure. The shear stress at the wall is shown in figure 
6. As seen from this figure, the flow starts from laminar 
and becomes turbulent at far downstream. The agreement 
between calculations and experimental data is reasonably 
good and the difference is less than 20% at x = 2.159 m. 
This difference seems to be due to the transition model 
and transition point. A better transition model, such as a 
precursor transition model, might improve the numerical 
solutions. 

The velocity and temperature profiles at x = 2.159 m 
are given in Figures 7 and 8. At this station, the flow is 
fully turbulent and the leading edge shock wave is located 
at y = 0.15 m. As seen from these figures, the calculated 
profiles agree relatively well with the test data. 

2) Hypersonic nozzle flow 

The Mach number contour for the hypersonic nozzle 
is given in Figure 9. The flow expands from sonic at the 
throat to about Mach 20 at x = .508 m as seen from this 
figure. The velocity and temperature profiles at x = .508 
m are shown in Figures 10 and 11. As seen from these 
figures, the profiles from present numerical solutions agree 
well with the data. The computed shear stress at the wall 
at x = 0.508m was 10.2l7V/m s , and the experimental data 
give the value of 10.60i\T/m 3 


Summary and Conclusions 

By solving the boundary layer equations for a flat 
plate at hypersonic speed, the variations of C cp and CkUb 
in the Baldwin-Lomax model with Mach number and wall 
temperature ratio were determined. By matching the solu- 
tions from Baldwin-Lomax model with those from Cebeci- 
Smith model, the new variations of the coefficients were 
obtained. The results show that the values of C cp and 
CkUb are functions of both Mach number and wall tem- 
perature ratio. C cp increases with Mach number and C cp 
for the adiabatic wall case is higher than that for the cold 
wall case. The variation of C*j e t shows that CkUb decreases 
with Mach number up to 5 and stays nearly constant when 
Mach number is greater than 5. The variations of C ep and 
CkUb were used for the calculations of hypersonic flows for 
both a 4° wedge and a nozzle. The solutions obtained by 
solving Navier-Stokes equations compared well with exper- 
imental data. 
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Figure 1 . - C cp variation with Mach number for a flat plate 
(adiabatic wall and R tz = 9.4 x 10 7 ). 



Figure 3. - C cp variations with Mach number and wall 
temperature ratio (adiabatic wall and ^ = 0.6, Tt = 
500° K, P t = 270. atm). 



MACH NUMBER MACH NUMBER 


Figure 2. - Cki,b variation with Mach number for a flat Figure 4. - Ckicb variations with Mach number and wall 
plate (adiabatic wall and R, x = 9.4 x 10 7 ), temperature ratio (adiabatic wall and ^ = 0.6, Tt = 

500’/:, P t = 270. atm). 




-COORDINATE (H) 


MACH 


NUMBER 


A 0.00 

B 2.13 

C *.27 

0 6.*0 

t 8 . 6 * 

F 10.67 

0 12.81 

H J *» 9 * 
1 17.00 

J 19.21 



Figure 9. - Mach number contour for a Mach 20 nozzle. 


CD 



Figure 10. - Velocity profile for a Mach 20 nozzle at x = 
0.508m. 


CO 



Figure 11. - Temperature profile for a Mach 20 nozzle at 
x = 0.508m. 


Table 1. Coefficients in Baldwin-Lomax Model 


Investigators 

c* 

Cud 

flowfield 

Baldwin and Lomax 1 

1.6 

0.3 

M=l, constant pressure 

York and Knight 5 

1.16-1.45 

.53-.65 

0 < M < 5, mass addition 

Granville 6 

1.25 

.59 

M « 1, constant pressure 

Granville 6 

.85-1.82 

.44-. 64 

M « 1, pressure gradient 
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